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ABSTRACT 

We measure radial gradients of the Mg 2 index in 15 E — E/SO and 14 5*0 
galaxies. Our homogeneous data set covers a large range of internal stellar ve- 
locity dispersions (2.0 < log o < 2.5) and Mg 2 gradients (AMg 2 / Alog r/r* e up 
to -0.14 mag dex" 1 ). We find for the first time, a noticeable lower boundary in 
the relation between Mg 2 gradient and a along the full range of a, which may be 
populated by galaxies predominantly formed by monolithic collapse. At high a, 
galaxies showing flatter gradients could represent objects which suffered either 
important merging episodes or later gas accretion. These processes contribute to 
the flattening of the metallicity gradients and their increasing importance could 
define the distribution of the objects above the boundary expected by the "clas- 
sical" monolithic process. This result is in marked contrast with previous works 
which found a correlation between AMg 2 / Alog r /r* and a confined to the low 
mass galaxies, suggesting that only galaxies below some limiting a were formed 
by collapse whereas the massive ones by mergers. We show observational evi- 
dence that a hybrid scenario could arise also among massive galaxies. Finally, we 
estimated A[Z/H] from Mg 2 and Hf3 measurements and single stellar population 
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models. The conclusions remain the same, indicating that the results cannot be 
ascribed to age effects on Mg 2 . 

Subject headings: galaxies: stellar populations — galaxies: formation — galaxies: 
elliptical and lenticular 

1. Introduction 

Two classical galaxy formation scenarios have been proposed to explain the build-up of 
elliptical galaxies: the Monolithic Dissipative Collapse (MDC) and the Hierarchical Clus- 
tering (HC). The existence of metallicity gradients in galaxies provide a key clue to help 
deciding between these two scenarios. 

In the MDC (e.g. Eggen, Lynden-Bell & Sandage 1962; Larson 1974; Chiosi & Car- 
raro 2002) a galaxy forms by means of a rapid gravitational collapse, with a considerable 
dissipation of energy from a cloud of primordial gas. During this process, star formation 
occurs in a short period of time. If the galaxy gravitational potential is strong enough to 
retain the gas ejected by stars and supernovae, it migrates to the galaxy's central regions. 
New generations of stars will be more metal rich than those in external parts and a steep 
radial negative metallicity gradient (MG) is established as a mass-dependent parameter. In 
this scenario ellipticals form at high redshift and on shorter timescales than spirals, and are 
assembled out of gas and not of pre-existing stars. Observational evidence for high mass 
galaxies already in place at higher redshifts have increased in the last years (e.g., Cimatti 
et al. 2004; Glazebrook et al. 2004). The appealing aspect of such models for the forma- 
tion of ellipticals is that they can explain most of the observational constraints relative to 
stellar populations. In particular, those assuming the inverse wind hypothesis (Matteucci 
1994) can also explain the observed increase of a/Fe with galaxy mass (e.g., Matteucci 2003; 
Thomas, Maraston, & Bender 2002) (see also Section 3). Observed MGs are reported in the 
literature by several authors (e.g., Gorgas, Efstathiou & Aragon-Salamanca 1990; Davidge 
1992; Carollo, Danziger & Buson 1993; Mehlert et al. 2003). Classical MDC models (e.g., 
Larson 1974) predict d(log Z)/d(log r) for ellipticals in the range —0.5 to —1.0, while the 
MGs correlate with global properties of elliptical galaxies, such that more massive galaxies 
have steeper MGs. The observational results of correlations between MGs and other internal 
global properties are controversial, as summarized by Kobayashi & Arimoto (1999). 

The scenario of HC suggests that present day galaxies form instead by a sequence of 
mergers of smaller objects (e.g., Steinmetz & Navarro 2002), and this is a natural consequence 
of the cold dark matter theory of cosmological structure formation. These models have 
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the advantage of working in a cosmological context, explaining the features of large scale 
structure. An important prediction of this kind of model is that the MGs should be erased by 
the merger process (e.g., White 1980). Flat MGs for some early-type galaxies were reported 
by Carollo, Danziger & Buson (1993). Additional observational evidence in favor of this 
scenario comes from the morphological and kinematical disturbances found in some ellipticals 
such as ripples/shells and multiple cores and kinematically decoupled cores (e.g.; Schweizer 
& Seitzer 1992; de Zeeuw et al. 2002). The morphology-density relation (Dressier 1980) can 
also be understood as the enhanced formation of ellipticals in high density environments, 
places where the efficiency of mergers is higher. In addition, the observed growth of the 
number density of red galaxies from z ~ 1.4 to the present time (Bell et al. 2004; Faber et 
al. 2005) argues in favor of the HC scenario. 

Kormendy (1989) considered a hybrid formation scenario for early-types where both 
pictures should be combined to explain observed structural properties of these types of 
galaxies. Recently, Kobayashi (2004) proposed the same idea based on simulations involving 
formation and chemodynamical evolution of galaxies. She shows that in a Cold Dark Matter 
scenario, the observed range of MGs requires the occurrence of both mechanisms described 
above. Kobayashi (2004) suggests that by using the MG values, it is possible to infer the 
merging histories of present-day galaxies. 

The goal of this letter is to investigate the dependence of MGs on a since this represents 
a powerful constraint on the different scenarios outlined above and is tightly linked to the 
mechanism of the formation of a galaxy. In Section 2, we present the sample of early-types 
used and a brief description about the method to measure MGs. In Section 3 we present 
and discuss our results, summarizing our conclusions in Section 4. 

2. The Sample and Line Index Determination 

Galaxies were selected from the ENEAR survey (da Costa et al. 2000), which contains a 
database of photometric (Alonso et al. 2003) and spectroscopic (Wegner et al. 2003) param- 
eters for a sample of early-type galaxies, which are representative of the nearby Universe. 
Several global parameters are available in the database such as magnitudes (m#), effective 
radii (r e ), mean surface brightness within r e , characteristic diameter (Dn), velocity disper- 
sion (a) and Mg 2 line indices in the Lick system. 

The spectroscopic data used in this letter were obtained with the 1.52 m telescope at 
the European Southern Observatory, during several runs. The sample consists of 15 E and 
E/SOs and 14 SO galaxies for which 1-D extractions were made as follows: the first aperture 
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(the central one) has 3 pixels (~2.5"). Successive lateral apertures are set in such a way that 
their central pixels are the outermost pixels of the previous apertures, keeping the same size. 
This process continues while a S/N > 20 is obtained for the spectra. When this condition 
fails, adjacent pixels towards the external parts of the galaxy are added iteratively, up to a 
maximum number of 5. The extraction process stops when it is not possible to keep S/N 
above 20, and this corresponds, typically, to a region of the object of ~ 0.5r e . No bias is 
expected in adopting variable radial zones as the gradients are essentially linear over the 
considered radial range. The final spectral resolution is ~ 4 — 5 A. The Mg 2 measurements 
were converted into the Lick/IDS system and brought to zero velocity dispersion following 
Worthey & Ottaviani (1997). 

A linear fit weighted by the index errors was used to measure its variation as a function 
of log(r/r*) : where r* is r e corrected for the galaxy ellipticity (r* = r e (l — e) -1 / 2 ). Because 
gradients are systematically flattened within the inner ~ 2" due to seeing effects, we ex- 
cluded these inner regions from the fits. We also discarded index measurements deviating 
by more than 2a from the linear fits. A detailed description of our data-reduction, analysis, 
and a comparison with gradient determinations by other authors will be the subject of a 
forthcoming paper. Meanwhile, a comparison of Mg 2 line gradients obtained in this work 
for 5 galaxies (IC 2035, NGC 2663, NGC 3557, NGC 3706 and NGC 5018) in common with 
Carollo, Danziger & Buson (1993); Carollo & Danziger (1994b) shows good agreement. The 
mean difference between these measurements is —0.005 ± 0.009 mag dex -1 which is quite 
reassuring. 

To perform the analysis described in the next section, we also included the data ob- 
tained by Carollo, Danziger & Buson (1993); Carollo & Danziger (1994a,b). To improve the 
homogeneity of the sample, we inspected visually the morphological classification for all the 
objects. Some of the Carollo, Danziger & Buson (1993) E and E/SOs were reclassified as 
SOs, based on visual inspection of DSS B,R and I available images. We indicate this when- 
ever necessary in the text and figures. As a final remark about the sample characteristics, 
we report that by means of a Kolmogorov-Smirnov test, the set of galaxies used in this letter 
has similar a distribution to that of ENEAR to a 99% confidence level. 

3. Results and Discussion 

To test one of the most important predictions of the MDC model - the dependence 
of MG with galaxy mass, we assume now that Mg 2 line index predominantly indicates 
metallicity. This line index cannot, in general, be taken as indicator of the metallicity only, 
since age effects should be considered. In fact, Mehlert et al. (2003) conclude that the 
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presence of a MG and the absence of gradients in age and a/Fe ratio for early- type galaxies 
in the Coma cluster imply that the intrinsic Mg — a relation is driven by metallicity alone. 
Nevertheless, we know that for our sample, made up of galaxies from field and other clusters, 
some more broader age distribution could play a role. While we are aware of these caveats, 
we first discuss Mg 2 gradients assuming this index is a good metallicity indicator. In order 
to further support our results, we combine measurements of Mg2 and Hf3 to infer [Z/H] 
gradients by comparison with the Thomas, Maraston, & Bender (2003) models (see below). 

The essential result of this letter is shown in panel (a) of Fig. 1 where we plot AM g 2 / Alog r/r* 
versus log a. Since there are few Es in the low a range, we also included SO galaxies, as- 
suming that they have similar metallicity properties as the Es (Bernardi et al. 1998). This 
figure shows that, while massive galaxies span a wide range of Mg 2 gradients there are no 
low mass galaxies with steep (very negative) Mg 2 gradients. We point out that the lower 
envelope of the data distribution seems to follow a linear relation between AMg 2 / Alog r/r* e 
and logo", which indicates that the formation of the galaxies occupying that locus might have 
been dominated by a monolithic collapse, as proposed by Kobayashi (2004). In her model, 
galaxies predominantly formed by monolithic collapse define a correlation with a extending 
to the low mass domain as suggested by the plot. Since our sample may be considered as 
representative of early-type galaxies in general, the lack of objects with steep gradients in the 
low a region is not a selection artifact, but instead a real feature of the galaxy distribution 
in the AM g 2 / Alog r/r* versus log a diagram. 

In panel (b) of Fig.l we plot the Mg 2 gradient versus log Mass (Mass oc r e a 2 ) using 
o and r e from our database. The distances were determined from the D n — a relation 
(Bernardi et al. 2002), adopting H = 75 km s _1 . This plot is very similar to that in panel 
(a), even though the Mass estimate reflects the uncertainties of both, a and r e , coupled to 
the hypothesis that the objects obey the virial theorem. 

Using only Mg 2 to represent metallicity, we may neglect possible age effects. Therefore, 
to make an estimate of the [Z/H] gradients, we used Mg 2 and Hf3 measurements and single 
stellar population models of Thomas, Maraston, & Bender (2003) for [a/Fe] =+0.3. This 
value of a/Fe was found by Mehlert et al. (2003) to be representative of all spatial regions 
of their sample galaxies. Metallicities were estimated for r e /8 and r e /2 from the Mg 2 and 
Hf3 values interpolated in the linear index vs. log r/r e relations. Galaxies for which one or 
both points fell out of the domain of the stellar population model predictions in index-index 
plots were excluded from this exercise. 

In Fig. 2 we display the A[Z/H}/ Alog r/r* versus a (panel a) and A[Z / H}/ Alog r/r* 
versus Mass (panel b) for 18 galaxies. We see that even though the position of some galaxies 
relative to the others changes substantially when we switch from Mg 2 to [Z/H], the overall 
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distribution of the data points in Fig. 2 is very similar to that in Fig.l. In particular, low 
mass galaxies appear to be characterized by low Z-gradients, whereas for high-mass galaxies 
there is a large spread in this quantity. As suggested by the simulations of Kobayashi (2004) 
there seems to be a lower envelope in the A[Z/H] versus Mass plane which coincides with 
the theoretical expectation from monolithic collapse models while the upper envelope for 
massive galaxies coincides with products of mergers. 

The fact that there are three data points with small positive gradients should not be 
taken at face value since they are consistent with a flat profile if one takes into account the 
associated errors. This should include errors such as H/3 emission contamination, variations 
on the true a/Fe ratios as well as the model uncertainties which are difficult to estimate. 
Therefore, we consider our error bars as an underestimate of the true error. Moreover, even 
if contrived, there are models that could explain positive gradients, at least for small mass 
objects (Mori et al. 1997). 

Gorgas, Efstathiou & Aragon-Salamanca (1990) and Davidge (1992), examining distinct 
samples of early-types, found a loose correlation between MG and a. They claim that the 
strong scatter in that relation is a consequence of different star formation histories. On the 
basis of a larger sample, Carollo, Danziger & Buson (1993) found that only low mass galaxies 
follow a MG — a correlation, while Kobayashi & Arimoto (1999) found no such correlation. 

Our results, on the other hand, show evidence of a dependence of the MG on a for the 
full range of a considered. In other words, the plot of MG vs. a now includes new galaxies 
for which a strong MG was found, providing a better definition of the lower envelope of the 
distribution for large a. It is interesting to note, that although Kobayashi (2004) simulations 
are within the hierarchical formation scenario, her results show that building galaxies from 
mergers of smaller objects with a high fraction of gas mimics a monolithic collapse. One 
difference of her simulations to a classic monolithic collapse is that the timescales for the 
formation of her monolithic cases are a bit longer. 

A quantitative comparison between the MGs measured here and those predicted by 
theoretical models (e.g., Kobayashi 2004; Pipino & Matteucci 2004) is made difficult by the 
fact that data and models sample different spatial regions. For instance, Kobayashi (2004) 
provide MG out to 2r e for objects in her simulations whereas our measurements extend only 
to ~ r e /2. Linear extrapolations of our measurements to larger radial distances cannot be 
trusted, as real galaxies could present substantial departures from linearity that would lead 
to large errors. Nevertheless, a qualitative comparison can be made between our results 
shown in Fig. 2, with those of simulations displayed in Fig. 6 of Kobayashi (2004). The lower 
boundary of the data points in our Fig. 2 seems to harbor the objects in Kobayashi's plot 
that were formed by monolithic collapse and minor mergers, while the objects distributed in 
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the upper region (weak MGs) may represent those with important merger histories. 

Pipino & Matteucci (2004) predict no correlation of MG with mass. In these models 
the MG developed in a galaxy depends on the interplay between star formation and stellar 
winds in the different galaxy zones. The flat dependence of the MGs with mass in this case 
can be explained by means of a strong "inverse wind effect" (Matteucci 1994). These models 
reproduce well the increase of the a/Fe ratio as a function of the galaxy mass by means of 
the inverse-wind hypothesis. Matteucci (1994) showed that to explain the increase of a/Fe 
in more massive galaxies it is important to assume that they complete their formation before 
the less massive ones. If this effect is strong enough, larger galaxies stop their star formation 
very fast preventing the development of strong gradients, leading to a flat dependence of MG 
with galaxy Mass (Pipino & Matteucci 2004). Another way in which models could produce 
both a-enhancement and stronger MGs in massive galaxies is via a variable initial mass 
function, whereby more massive galaxies would have IMFs biased towards high mass stars 
(Matteucci 1994). 



4. Summary 

We analyzed MGs for a sample of early-type galaxies and our results show: (a) at 
least part of the objects follow the predicted MG — a relation claimed to exist in the MDC 
scenario, defining a lower boundary in a MG vs. a diagram; (b) the remaining objects may 
represent a sequence of increasing relative importance of the HC process. Particularly, higher 
mass galaxies scatter vertically in this diagram, indicating a contribution by both formation 
processes. These results give observational support to the hybrid scenario proposed by 
Kormendy (1989) and Kobayashi (2004) and show that we can interpret galaxy formation 
not as an exclusive dominance of either MDC or HC scenarios, but instead, both mechanisms 
contribute to the origin of these objects. Each galaxy has its particular formation history 
depending on the merger events of its building blocks, their nature (predominantly gas or 
stars) and their efficiency to collapse. In this case the location of a galaxy in a MG — a 
diagram may be a useful tool to infer the relative importance of mergers or collapse to its 
formation. 

Additional observations underway will be used to increase our statistics and to cover 
more homogeneously the range of parameters examined. This will allow us to investigate 
the existence of a correlation between MGs and environmental density and to examine the 
radial behavior of a/Fe, to be presented in a future paper. 
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Fig. 1. — Radial Mg 2 gradients as a function of central velocity dispersion a (panel a) and 
versus Mass (panel b) for E, E/S0 and 5*0 galaxies. We also include the data for the 
early-types compiled by Carollo, Danziger & Buson (1993); Carollo & Danziger (1994a,b) 
(CDB93+CD94ab). A few Es of these authors turned out to be SOs. Their Es and E/SOs 
are indicated by "V" and the SOs by "A". The symbol "V" in the label inside the figure for 
CDB93+CD94a,b E, E/SO contains an average error bar for their data. The 5 galaxies in 
common with their compilation are connected by a solid line to give an idea of the differences 
between our and their determinations of AMg 2 and a. The size of the symbol is proportional 
to the fraction of r e in which we measured the AMg 2 (the bigger ones means r/r e » 1). 
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Fig. 2. — Radial A[Z/H] as a function of galaxy central velocity dispersion a (panel a) 
and Mass (panel b). The number of points is smaller in this plot because some line indices 
measurements of Mg 2 and HP used to determine [Z/H] fall outside the grid of single stellar 
population models by Thomas, Maraston, & Bender (2003). 



